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Abstract 

We have previously proposed a light pseudoscalar Higgs boson in the next-to-minimal supersymmetric 
standard model (NMSSM), the A^, as a candidate to explain the HyperCP observations in pfi~^i_i^ . 

In this paper we calculate the rates for several other rare decay modes that can help confirm or refute this 
hypothesis. The first modes we evaluate are — > ttttA^, -which are interesting because they are under 
study by the KTeV Collaboration. We next turn to — > ttttA^, -which are interesting because they are 
independent of the details of the flavor-changing sector of the NMSSM and may be accessible at DA^NE. 
For completeness, -we also evaluate il~ — > H~j4'j*. 
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I. INTRODUCTION 



The HyperCP Collaboration recently observed three events for the rare decay mode S'^ pjj,~^iJ,~ 
with dimuon invariant masses narrowly clustered around 214.3 MeV [1]. It is possible to account 
for these events within the standard model (SM) when long-distance contributions are properly 
included [2, 3]. However, the probability that the three events have the same dimuon mass, given 
the SM predictions, is less than one percent. This result has prompted several studies investigating 
the consequences of a new state with this mass [4-6] . 

In particular, it was pointed out that the flavor-changing coupling of the new state to ds has 
to be (dominantly) of a pseudoscalar or axial-vector nature to explain why it has not been seen 
in K ^ 7r/x"'"/x~. This would still allow the new particle to be observed in the other rare modes 
K nn fi~^fi~ and Q~ S'fi'^fi''. Predictions for the new particle's contributing to these modes, 
consistent with existing constraints, were made in Refs. [4, 5]. These predictions indicate that there 
could be evidence for the particle in the data already taken by the KTeV Collaboration, speciflcally 
in the mode Kl — > Tr^Tr^//"*"//" currently being studied [7]. 

Beyond the above-mentioned theoretical analyses, to explore the possible consequences of the 
HyperCP result in greater detail one has to incorporate some model dependence. To this end, 
various ideas have been proposed in the literature [8, 9]. Specifically, we have demonstrated that 
a light pseudoscalar Higgs boson in the next-to- minimal supersymmetric standard model (NMSSM), 
the A^, could be identified as the possible new particle responsible for the HyperCP events while 
satisfying all constraints from kaon and E-meson decays [9] . 

Now, it is long known that kaon decays involving a light Higgs boson, such as in the NMSSM, 
receive two types of contributions that can be of comparable size: two-quark contributions in which 
the flavor change occurs in one- loop processes involving the light ^4^, and four-quark contributions 

in which the flavor change occurs via a tree-level standard-model W exchange with the light 
radiated off one of the light quarks [10, 11]. Not too long ago, we showed that the same situation 
occurs in the case of light Higgs production in hyperon decays [12]. 

In this paper we revisit the modes K nnAi and Q~ 'E~A^ in order to present a complete 
prediction within the model suggested in Ref. [9], plus the possible modifications recently pointed 
out in Ref. [13]. This differs from the model- independent studies of Refs. [4, 5] in two important 
ways. Within the NMSSM, we can identify the effective scalar and pseudoscalar couphngs of the 
model-independent studies with specific one-loop processes. Here we consider not only the chargino- 
mediated diagrams of Refs. [9, 14], but also the gluino- and neutralino-mediated diagrams discussed 
in Ref. [13]. In addition, we include the four-quark contributions which are missing in Refs. [4, 5]. 
These four-quark contributions were shown in Ref. [9] to be essential to evade the bounds arising 
from the nonobservation of the A^ in K ^ nn^n^ modes [15]. Following our earlier work [4, 9], 
we will here assume that B{A'( — > jjf^jJT) ~ 100%, 

Additional processes where such a light A^ would appear have been recently studied in the 
hterature: collider signatures for a light A\ [16], S-meson decays [17], and radiative quarkonium 
decays [18] . The latter are especially useful because, being flavor conserving, they are independent 
of the specifics of the one-loop fiavor-changing couplings and follow directly from the tree-level 
couplings of the A^ to down-type quarks. 
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Finally, in this paper we also consider the modes 77 — > ttttAi which, like radiative quarkonium 
decays, are flavor diagonal and only sensitive to the tree-level couplings of the A^. The prediction 
for these modes is, therefore, much less model-dependent. We find a rate two orders of magnitude 
larger than the corresponding SM rate, that could be probed at DA$NE. 

II. THE LIGHT ^5 IN THE NMSSM 

In this section, we briefly review some features of the NMSSM that are relevant to our study. 
The model is an extension of the minimal supersymmetric standard model (MSSM) and provides 
a solution to the so-called /x-problem of the MSSM [19]. In the NMSSM, there is a gauge-singlet 
Higgs field N in addition to the two Higgs fields and responsible for the up- and down-type 
quark masses in the MSSM. As a result, the physical spectrum of the extended model has two 
additional neutral Higgs bosons: one a scalar and the other a pseudoscalar. 

We follow the specific model described in Ref. [14], with suitable modifications. The superpo- 
tential of the model is given by 

W = QY^H^U + QY,H,D + LY^H.E + XH.H^N - ^kN' , (1) 

where Q, U, D, L, and E represent the usual quark and lepton fields, l^„,rf,e are the Yukawa couplings, 
and A and k are dimensionless parameters. The soft-supersymmetry-breaking term in the Higgs 
potential is 

Koft = rnljHj' + mljHf + m%\N\' - {XA^H.H^N + \kA,N' + H.c.) , (2) 

and the resulting Higgs potential has a global U{1)r symmetry in the limit that the parameters 
A^,Ak^Q [20]. 

The NMSSM has two physical CP-odd Higgs bosons which are finear combinations of the pseu- 
doscalar components in Hd, and N in the model mix, with the ^ being the lighter mass- 
eigenstate with mass given by 

m\ = ?,kxAj. + C(l/tan/?) (3) 

in the large-tan /3 limit, where x = (N) is the vacuum expectation value of A*" and tan /3 is the ratio 
of vacuum expectation values (VEVs) of the two Higgs doublets. If the U{1)r symmetry is broken 
slightly, the mass of A^ becomes naturally small, with values as low as ~100 MeV phenomenologi- 
cally allowed [14, 20, 21]. 

In the large-tan /9 limit, the A^ is mostly the singlet pseudoscalar and couples to SM fields 
through mixing. Also in the large-tan /3 limit, its couplings to fermions are suppressed by a factor 
of tan P with respect to those of the A'^ in the MSSM [14, 20] . In particular, this makes the tree-level 
couplings to up-type quarks negligible. The tree-level couplings to down-type quarks and charged 
leptons can be described in terms of one parameter. 
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where the parameter l^, involves both the different Higgs VEVs and soft-supersymmetry-breaking 
parameters, 



V 5. 



(5) 



with V — 246 GeV being the electroweak scale and S_ = {Ax — 2kx)/{Ax + kx). Requiring the 
mass of the heavier pseudoscalar not to exceed 500 GeV, Ref. [14] finds a lower bound ^0.1 
for tan /3 = 30. At the same time, the contribution of to the muon anomalous magnetic moment 
results in an upper bound ^ 1.2 [4]. 

As shown in Rcf. [9], this scenario leads to four-quark contributions that easily reproduce the 
HypcrCP result. Unfortunately, they are also in conflict with the nonobservation of the in 
K nfi^/j^ [15]. To satisfy these bounds, it is necessary to include contributions from one- 
loop flavor-changing diagrams in the NMSSM, which depend in a complicated way on the many 
parameters of the model. We can describe them in terms of an effective Lagrangian for the sdAi 
couplings as 



In Ref. [9], we followed Ref. [14] to consider only chargino-mcdiated one-loop diagrams in the 
large-tan /3 limit. Furthermore, we selected the supersymmetric parameters so as to suppress the 
6 — > s transition and in this way satisfy the nonobservation of A^ in B decay [22] . This scenario led 
to Cj^ — —Cj^m^/nig — —2gj^mjv with Qj^^ ~ 10~^ [9]. More recently, Ref. [13] has pointed out 
a different scenario in which Cl,r also receive contributions from gluino- and neutralino- mediated 
one-loop diagrams. Although the gluino-mediated contributions are suppressed by a factor of tan 13 
compared to the chargino contributions, the former are proportional to the strong coupling 
compensating for the suppression factor, and hence can be as important as the latter. Moreover, 
in some regions of the parameter space the neutralino-mediated contributions could be comparable 
to the gluino-mediated ones [13]. If all the different contributions are similar in size, then Cl,r can 
become effectively independent. 

This opens up the possibility of satisfying the kaon bounds without the four-quark contributions 

by having Cl ^ —Cr, which results in an effective sdA^ coupling that is mostly pseudoscalar. The 
HyperCP observation can then be explained as in the model- independent analysis of Refs. [4, 5]. 
However, to have Cl ^ —Cr requires some sort of fine tuning. Furthermore, the four-quark 
contributions may not necessarily be negligible. In our analysis, we will thus keep Cl and Cr 
independent and constrain them with data. Also, we will assume that CP is conserved and hence 
Cl,r are real. 



iC - iC 



(6) 



where the parameters Cj^ ^ are generally independent. 
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III. |A5| = 1 DECAYS 



A. Two-quark contributions 



To evaluate hadronic amplitudes induced by the sdA^ interactions, we employ chiral perturba- 
tion theory (xPT). Thus, the leading-order chiral realization of Cj^sd above is £_4 in Eq. (Al) in 
Appendix A, which also contains other relevant chiral Lagrangians. From £_4 and the chiral strong 
Lagrangian Cs in Eq. (A3), we derive the leading-order diagrams shown in Fig. 1 for K ttttA^ 
and S^A^. The resulting amplitudes are 



TC 7C 



B, {C, - C^) m\_^ - ml 



m 



A 



(7) 



^0 {Cl - Cr) m\-m\- m^o^o 



4^2/ 



m\-m\ 



(8) 
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where m|-y = {px +PyY- The same Lagrangians also yield 



TX 



[CI + Q) 



(9) 
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2 m]^ 
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previously derived in Refs. [12, 13]. Hence we also adopt D — F = 0.25. 



(11) 
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9 ^1 



K 



n ^0 



n- 



(b) 



(a) 

FIG. 1: Diagrams contributing to (a) K mrA^ and (b) — 
in xPT. The square vertices come from in Eq. (Al), and the solid dots from £g in Eq. (A3). 



'A^ arising from -C^^^^ at leading order 
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B. Four-quark contributions 



From given in Appendix A, we obtain the leading-order diagrams shown in Figs. 2 and 3 
for K irnA^ and Q~ S~y4°. The resulting amplitudes are 

8 

M,^{K'^n^n-A[) = J^^t , (12) 

i=l 

A^,,(K°-7rVA?) = (13) 

i=l 

^ u 

where the expressions for A4^^'^^ and fe7r,r),j7' have been collected in Appendix B, and B^-^^o is 
related in xPT to the dominant P-wave amplitude for S^vr*^ by S^tt^) = 

iB-^-T^o ij)^) ^u^u'^. Hence the Vt^ H^vr^ data yields B^-^^o = —8.17 x 10^''. We note that 
the 78 contributions to M.^^[K n-nA^ cancel completely, which is expected due to the fact that 

the 7g terms in could be rotated away if the baryonic part were absent [11]. We also note that 
the 7g contribution to Vt~ S"??*^ appears only at next-to- leading order. 

For K — > T[A\ and S+ — > pA\, the four-quark amplitudes were previously calculated in Ref. [12]. 
For = 0, they can be rewritten as^ 

- {b^sg - b^,cg) [m]^ - ml) 7* , (15) 
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FIG. 2: Diagrams contributing to K ^ ttttA^ arising from four-quark operators, where V = TT^,r],r]'. 
The dots come from ci"^^ in Eqs. (A3) and (A5), whereas the square vertices are from in Eqs. (A4) 
and (A6). 



^ There is a typo in the last line of Eq. (70) in Ref. [12]. The term —{l^ + l^)m^ should be corrected to — (3Z^ + ^„)mJ 
This error, however, did not occur in our computation. 
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Q- — 

FIG. 3: Diagram contributing to E^A^ arising from four-quark operators. 

M,^{K'^7c'A[) = '-^\3b^ {2ml - ml - m'^) " [b^c^ + b^,s,) {2ml + ml - Sm^,) 



+ VS {b^s, - b^,c,) {ml - ml)h; I, , (16) 



^pA) = ^ {-b^ + b^c, + b,,,s,) tp {A^^o - 5,^075) S+ (17) 

where Ap^o = -3.25 x 10"'^ and Bp^o = 26.67 x 10"^ up to an overall sign, extracted from 
S+ pTT^ data. 



C. Total contributions 

In this section we present numerical results for the different modes including all contributions 
to the respective amplitudes. We begin by determining the region in the {Cl+Cr, Id) parameter 
space that is allowed by both the n'^fi'^fi' and Ks — > 7r°/i+/i~ constraints. We show this 

in Fig. 4. Notice that only small values of C^+C/j are allowed. This corresponds to the conclusion 
of the analyses of Refs. [4, 5] that the effective sdA^ scalar coupling is severely constrained by these 
decay modes. That case, without the four-quark contributions, corresponds to Id = in this plot. 
The inclusion of the four-quark contributions does not change this conclusion, but simply shifts the 
allowed region due to the interplay between the two- and four-quark contributions. 

For definiteness, we select = 0.35 as in Ref. [9] and study the allowed region in the 
{Cl+Cr, Cl—Cr) parameter space. We display in Fig. 5 the lightly shaded (yellow) region that 
reproduces the HyperCP result for S"*" — >■ pfi^fi" (at the one-sigma level combining statistical and 
systematic errors in quadrature). The darkly shaded (red) vertical band covers the region that 
satisfies the constraints from the nonobservation of A^ in K ^ 7r/i+yU^ modes. For = 0, these 
(yellow and red) areas would both be centered at the origin. The (black) overlap between these 
regions is the allowed parameter space that we use for our predictions. Also displayed on the ver- 
tical band is an unshaded (white) thin area corresponding to the = —Cj^m^/m^ = —2g_^m^/v 
scenario of Ref. [9]. 

With these results, we show in Fig. 6 the predicted branching ratios (solid curves) for 
TT+TT-AO and tt^ttM? as functions of Q-C^ for = 0.35 and Q + = 4 x 10"". The 

range of each of these predictions over the allowed values of Cj^—Cj^ is larger than that obtained 
in Ref. [4] , due partly to the presence of the four-quark contributions and partly to the uncertainty 
in the HyperCP measurement. Each of the solid curves has a minimum that is not zero, as the 
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FIG. 4: Regions in the {Cl+Cr,Iii) parameter space allowed by vr^/U^/i^ (shaded, blue) and 

Kg — > TT^fi'^fi^ (lightly shaded, green). The overlap (dark, red) band covers points that satisfy both 
constraints. 
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FIG. 5: Regions in the {Cl+Cr^Cl—Cr) parameter space reproducing the HyperCP result for S+ — > 
pfi'^fi^ (lightly shaded, yellow) and respecting the K — > ■Kfi'^fi^ bounds (darkly shaded, red) for 1^ = 0.35. 
The overlap (black) areas cover points satisfying both the hyperon and kaon constraints. The unshaded 
(white) region on the vertical band corresponds to the special case discussed in Ref. [9]. 



two- and four-quark contributions have different kinematical dependences and hence do not cancel 
in general. The rates for most of the allowed regions are significantly large, but those around the 
minima may be too small to be observed. For comparison, we also show in Fig. 5 dotted curves 
representing the branching ratios obtained from the two-quark contributions alone and vertical 
(green) dashed lines each indicating the narrow range of C^— found in the scenario of Ref. [9]. 
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FIG. 6: Predicted branching ratios (solid curves) for Kj^ ir^ir^Ai and tt^tt^A^ with = 0.35 

as functions of C^—C^. The dotted curves result from the two-quark contributions alone. The shaded 
(pink) bands indicate the allowed ranges of C^—Cj^ as determined from Fig. 5. Each vertical (green) 
dashed line corresponds to the special case discussed in Ref. [9]. 

Since the values of Ap^^o and Bp^o in Eq. (17) are determined from experiment only up to an 
overall sign, we should also consider the possibility that the two and four-quark contributions to 
—> pii'^n~ have a different relative sign. This yields a different allowed range of Cl—Cr, as 
can be seen in Fig. 7. We display the resulting predictions for K^^ ttttA^ in Fig. 8, whose ranges 
over the allowed regions turn out to be roughly only half as large as those in Fig. 6, respectively. 

Finally, from the results of Fig. 5 we display the predicted branching ratio (solid curve) for 
E^A^ in Fig. 9. The range of the prediction over the allowed values of Cl—Cr is again 
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FIG. 7: The same as Fig. 5, except that the relative sign between the two- and four-quark contributions 
to S+ — > pij,~^ii~ is the opposite. 
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FIG. 8: The same as Fig. 6, except that the allowed range of Cj^—C^ is from Fig. 7. 

larger than that obtained in Ref. [4] due to the presence of the four-quark contributions as well 
as to the experimental error. The best limit for this mode currently available was reported by 
HyperCP [23], BiVL^ E^^^^^) < 6.1 x 10^^ at 90% C.L., whereas the standard-model prediction 
is Bsm{^^ A^^Ai ) = 6.6 x 10^^ [24]. Therefore, the Q~ —>■ S'A^ rate for most of the allowed 

regions is substantial, but the curve has a zero, around which the rate is too small to be observed. 
The significant enhancement possible with respect to the SM rate lends support to pursuing a future 
Q~ experiment [25]. For comparison, we also display in Fig. 9 the dotted curve representing the 
branching ratio obtained from only the two-quark contributions and the vertical (green) dashed 
lines corresponding to the special case of Ref. [9]. In Fig. 10 we show the corresponding prediction 
with the allowed range of Cl—Cr from Fig. 7. 




10'°(Q-C«) 



FIG. 9: Predicted branching ratio (solid curve) for Q.~ — > A\ with = 0.35 as function of Cj^—Cj^. 
The dotted curve results from the two-quark contributions alone. The shaded (pink) bands indicate the 
allowed ranges of C^—Cj^ as determined from Fig. 5. The vertical (green) dashed line corresponds to the 
special case of Ref. [9]. 
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FIG. 10: The same as Fig. 9, except that the allowed range of C^—Cj^ is from Fig. 7. 



IV. FLAVOR-CONSERVING DECAYS r? ttttA? 



These modes are special because they involve only flavor-diagonal interactions. As such, they 
are not sensitive to the unknown parameters in the flavor sector of the model that give rise to the 
two-quark amplitudes. The predicted rates follow only from the tree-level diagonal couplings of 
and in this way they are similar to the radiative quarkonium decays proposed in Ref. [18]. These 
7] decays are also analogous to the rj decay with a light CP-even Higgs boson which was severely 
constrained by data [26]. 

The leading-order amplitude for t] — ^ ttttA^ comes from the two diagrams in Fig. 11. It is the 
same for rj 7r"'"7r~y4° and f] n^n^A^, 



Miv^nirA',) = ^^^^3{cg-V2sg) + 6^ (l - c,s, + s^) + b^, {^2 - CsS^ - Vs si) I 



d 



For an rj-r]' mixing angle of = — 19.7°, we then find 

B{r] 7r+7r-A°) = 5.4 x 10"^/^ , (19) 
B{t] ^ ttVA?) = 3.2 X 10"'/^ . (20) 

Allowing the mixing angle to vary between —15° and —25° would result in 20% changes. The rate 
for the neutral-pion mode is not exactly half the rate for the charged-pion mode because we have 
used physical masses for the numerical estimate. 



FIG. 11: Diagrams contributing to r] ^ ttttA^ induced by flavor-diagonal couplings of to light quarks. 
The dots come from ci"^^ in Eqs. (A3) and (A5), 
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The best limit currently available on any of these modes comes from the CELCIUS/WASA 
collaboration. At the 90% C.L. they find [27]. 

B{r] 7r+7r-//V~) < 3.6 x 10"^ . (21) 

Presently this docs not place a stringent bound on the coupling la, giving \ld\ < 26. Nev- 
ertheless, Eq. (19) is a very significant enhancement over the expected standard-model rate, 
Bsuiv 7r+7r-//+/i-) = {T.Bttr) x 10"'' [28, 29], and may be accessible to DA$NE [30]. 

V. SUMMARY AND CONCLUSIONS 

We have studied several rare decay modes involving a hght CP-odd Higgs boson in the NMSSM. 
In the analysis, for the flavor-changing modes, we have consistently included the two-quark contri- 
butions in which the flavor change occurs in one-loop processes involving the light and the four- 
quark contributions in which the flavor change occurs via a tree-level standard-model W exchange 
with the light radiated off one of the light quarks. The interplay between these two contributions 
was crucial to evade the bounds arising from the nonobservation of the A^ in K ^ 7iii~^ii~ modes 
in our previous analysis [4]. 

For the two-quark contributions, we have considered a somewhat general scenario in which the 
coefficients Cl,r are effectively independent. We have started with the large-tan /3 hmit where 
chargino-mediated one-loop diagrams dominate, but we have also allowed for the possibility of 
having sizable neutralino- and gluino-mediated one-loop diagrams. In this more general scenario, it 
would also be possible to evade the K vr/x^/i" bounds even if the four-quark contributions were 
absent. 

We have evaluated the rare modes Kj^ titiAi which depend on both the two- and four-quark 
contributions. We have found that their rates are significant for most of the allowed parameter 
space. These modes are of immediate interest because they can be studied with KTeV data. It is 
expected that these studies can help confirm or refute the hght A^ hypothesis as a candidate to 
explain the HyperCP events in E''" — > p^^ ^~ . 

We have also studied the modes rj — > ttttA? which depend only on the tree-level couphngs of 

the A\. Therefore, the predictions for these modes are much less model-dependent and should 
be of interest for future experiments at DA$NE. In particular, the 745-mediated contribution to 
7] 71^ 7r~ fi~^ fi~ can be much larger than the SM contribution. 

Finally, we have revisited the mode Q~ S'A^ S'fi^fi^ to include both the two- and 
four-quark contributions. We have found that its decay rate could be substantially enhanced with 
respect to the — > S~/i+//~ rate in the SM. This should give additional motivation for conducting 
experimental studies on the in the future. 
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APPENDIX A: CHIRAL LAGRANGIANS FOR VARIOUS INTERACTIONS 

The Lagrangians we have collected here contain not only the baryon- and meson-octet fields, but 
also the baryon-decuplet fields. Since we already derived or used some of the following formulas in 
Refs. [4, 12], further details can be found therein. 

The chiral realization of C^sd in Eq. (6) can be obtained employing the operator matching of 
Ref. [4]. Thus at leading order 

£^ = hi,{B {hj„B}) + h^{B[hj,,B]) + b,{hj,) {BB) + \fB^{hj,) 

+ cf^hXr^ - Co {h^) f'^T^ + H.c. , (Al) 

where f = f^ = 92.4 MeV, B^ = 2031 MeV, and 

hA = -iiCR^he + C^iK) A\ . (A2) 

To derive amplitudes, we also need the chiral Lagrangian for the strong interactions of the 
hadrons [31, 32]. At lowest order in the derivative and expansions, it can be expressed as 

4 = {Bir{d,B+[V^,B]))-m^{BB) + D{B^^^,{A^,B]) + F{Bri,[A^,B]) 

+ 6^ {B {M+, B}) + hp {B [M+, S]> + 6o {M^) (BB) + \f {d^E^ d^E) + \fB^ (M+) 
- iTpT^ + rriT Ti'T^ + C {TI'A^B + BAJTi") + c T^M+T^ - (M+> T^T^ , (A3) 

where only the relevant terms are displayed and \C\ = 1.7. 

The leading-order Lagrangians relevant to the four-quark interactions involving the A^ can be 
derived, following the prescription described in Refs. [11, 12], from above and from the mass 
term in the leading-order Lagrangian for the |A7| = | component of the effective Hamiltonian 
transforming as (8l, Ij^), namely [33] 

= hj,{B{eH, B]) + hp{B [i'^hi, B]) + ^J^hd^T.d^T.'^) + 2%pB,{hiM^e) 
+ hcf'^eKT^ + H.c. , (A4) 

where 73 = —7.8 x 10~^. Thus we have 

= (bj,{B{M_,B}) + bp{B[M_,B]) + b,{M_){BB) + \fB^{M_) 

\ iA^ 

+ c™_7; - Co <M_> T'^T^j , (A5) 

40 

= '2%fB,{hCM_e)-l + H.c. , (A6) 
13 



where M_ = C^M^^ - ^M^^, with M = diag(0, l^rh, l^m^) at large tan/3. We include the SU(3) 
singlet rji in Ci"^ by replacing E with E exp(i-y/2/377^//) and adding the anomaly-generated term 



^^2 



(A7) 



which modifies the rj^-A'^ mixing generated by Cf. The physical rj and rj' fields are related to 77^ 
and the SU(3) octet 773 by rj — rj^ Cq — rji Sg and 77' — r]^Sg + rji Cq, where = cos 9 and — sin ^. 
As in Ref. [12], we adopt ttiq = 819 MeV and 9 = -19.7°. 

APPENDIX B: FOUR-QUARK CONTRIBUTIONS TO K'^ ttttA? 



In the case of ^ n'^n A^, for 1^ — the eight diagrams in Fig. 2 yield, respectively. 
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where m^y = (p^ +py)2, 
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In the case of ^ tt^tt^AJ, the diagrams in Fig. 2 yield, for = 0, 
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In the expressions for M.'l' or above, we have kept the terms proportional to 7g in order to 
check our algebra. As explained in Ref. [11], the 73 terms in can be rotated away for kaon 
decay, and we have verified that the 7g terms cancel accordingly in the sum of the contributions. 
In our numerical evaluation, 73 is thus set to zero. 
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